%ROTATION_MATRIX
%
%INPUT
%alpha - rotation angle
%
%OUTPUT
%rotation matrix
%
function matrix = rotation_matrix(alpha)
s = sin(alpha);
c = cos(alpha);
matrix = [c -s;s c];